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Abstract 

It  has  been  proposed  that  source  localization  in  passive 
sonar  be  accomplished  via  some  form  of  Matched  Field 
Processing  (MFP).  Full-coherent.  Minimum  Vr.riance 
Distortionless  Response  (MVDR)  MFP  assumes  complete 
and  perfect  characterization  of  the  underwater  multipath 
propagation  channel  and  is  known  to  be  extremely  sensi¬ 
tive  to  mismatch  between  model  parameters  and  actual 
environmental  parameters.  We  present  a  minimax. 
multi-rank  signal  variation  of  MVDR  MFP  referred  to  as 
Semi-coherent  MVDR  MFP.  Simulation  results  are 
presented  demonstrating  Semi-coherent  MVDR  MFP  to 
be  both  relatively  robust  to  mismatch,  with  respect  to 
relative  amplitudes  and  phases  amongst  multipath 
arrivals,  and  comparable  in  performance  to  Full- 
coherent  MVDR  MFP  under  no  mismatch  conditions. 

1.  Introduction 

Consider  a  candidate  source  location  (i.  e„  a  point 
on  a  search  grid)  designated  by  the  position  vector 
r*=  (R ,  0 ,  z)  with  respect  to  a  cylindrical  coordinate  sys¬ 
tem  centered  at  some  reference  point  in  the  array. 
Through  acoustic  ray  tracing,  we  determine  the  respec¬ 
tive  arrival  angle  of  each  ray  path  between  the  candidate 
source  location  and  some  reference  point  in  the  array. 
Let  L  denote  the  number  of  dominant  ray  paths  based  on 
ray  tracing  prediction.  In  the  case  of  a  linear  array,  we 
denote  the  Lxl  vector  composed  of  the  L  conical  arrival 
angles  as  0,  i.  e.,  0  =  (0j  ,  02  .  * '  ‘  .  0l).  where  0,  is  the 
conical  arrival  angle  associated  with  the  i-th  ray  path, 
i=l,...L-  With  further  modeling,  we  could  determine  the 
relative  amplitudes  and  phases  amongst  the  L  multipath 
arrivals.  However,  the  relative  phases  of  the  multipath 
arrivals  can  change  dramatically  with  small  changes  in 
the  ocean  parameters,  a  change  in  the  ocean  depth,  for 

*Thi»  work  w»l  lupporled  by  the  Nival  Ocean  Systems  Censer 
(NOSQ  under  contract  no.  N66001-87-D-0I36  with  Office  of 
Naval  Technology  (ONT)  funding. 
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example  [1].  As  a  step  towards  developing  a  robust  pro¬ 
cedure,  we  will  not  assume  knowledge  of  the  relative 
amplitudes  and  phases  of  the  multipath  arrivals  for  any 
source  location.  Rather,  we  will  only  assume  knowledge 
of  0,  and  the  corresponding  steering  vector,  denoted  a,, 
i=l,,„JL.  for  cadi  of  the  L  dominant  ray  paths. 

2.  Semi-coherent  MVDR  MFP 

Semi-coherent  MVDR  MFP  is  developed  as  a 
minimax  approach  to  the  source  localization  problem. 
Let  N  denote  the  number  of  sensors  comprising  the  array. 
An  Nxl  weight  vector  is  constructed  to  minimize  the 
average  power  of  the  corresponding  beam  former  output 
under  a  constraint  on  the  gain  and  phase  response  in  each 
of  the  L  multipath  arrival  directions.  The  gain  and  phase 
response  pairs,  one  for  each  multipath  arrival  direction, 
are  jointly  determined  so  as  to  maximize  the  SNR  of  the 
beamformer  output.  The  minimization  stage  may  be 
mathematically  posed  as  the  following  constrained 
optimization  problem: 

Minimize  w^r*  )R„wCr  )  (2.1) 

w(7) 

subject  to:  AH(0)w(7)  =  8 

where  5  =  [Sj  ,  82  ,  •  •  •  .  5l)t,  an  Lxl  vector,  and 
A(8)  =  [at  ,  82  >  •  •  •  ,  aL],  an  NxL  matrix.  The  magni¬ 
tude  and  phase  of  6;  represent  the  gain  and  phase 
response,  respectively,  in  the  i-th  multipath  arrival  direc¬ 
tion,  i=l,...,L.  Assuming  L  <  N,  AH(0)w(7)  =  8 
represents  an  underdetermined  system  of  equations.  The 
solution  to  (2.1)  is  simply  the  minimum  norm  solution  to 
AH(0)w(7)  =  5  in  a  Hilbert  space  with  inner  product 
<x  ,  y>  =  y^x: 

w(7)=  R„A(0)  |ah(0)R^1  A(0)j  ‘5  (2.2) 

Motivated  by  the  work  of  Krolik,  Lynch,  and  Swingler  in 
developing  "Incoherent"  MVDR  MFP  [l)1 ,  we  define  the 
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value  of  the  Semi-coherent  MVDR  MFP  ambiguity  sur¬ 
face  at  a  point  7  for  a  given  beam  response  constraint 
vector  5,  denoted  Ssemif?  ;  5),  to  be  the  SNR  of  the 
beamformer  output  obtained  with  w(7)  in  (2.2).  Actu¬ 
ally,  since  we  don’t  have  access  to  signal-only  data,  we 
will  define  it  to  be  the  signal  plus  noise  to  noise  ratio. 


For  the  sake  of  simplicity,  we  will  assume  the 
noise  to  be  spatially  white  with  the  power  of  the  noise  at 
each  sensor  element  equal  to  a2.  h1  this  case,  the  output 
signal  plus  noise  to  noise  ratio,  (S+N)/N,  is  proportional 
to  wH(7)Ruw(7)/  wH(7)w(7).  With  w(7)  given  by 
('-.2),  Ssani(?  5)  may  be  expressed,  after  some  mani¬ 
pulation,  as 


Ssenu<7:5) 


5hM7'M1M715 

5hM71M2M7'5 


(2.3) 


subject  to:  wH(7)  (A(8)c(7)1  =  1  (3.1) 

Thus,  in  addition  to  knowledge  of  the  arrival  angle  and 
corresponding  steering  vector  for  each  ray  path.  Full- 
coherent  MVDR  MFP  also  assumes  knowledge  of  the 
relative  amplitudes  and  phases  amongst  the  L  ray  paths. 
In  addition.  Full-coherent  MVDR  MFP  assumes  the  mul¬ 
tipath  arrivals  to  be  100%  correlated.  In  contrast,  proper 
operation  of  Semi-coherent  MVDR  MFP  is  not  premised 
on  such  --  the  multipath  arrivals  corresponding  to  the 
same  source  may  be  partially  correlated  or  100%  corre¬ 
lated. 

Similar  to  the  solution  to  (2.1),  the  solution  to  (3.1) 
is  the  minimum  norm  solution  to  the  constraint  equation 
with  the  norm  induced  by  the  inner  product 
<x  ,  y>  =  y“R„x. 


where  M,  and  M2  are  defined  as 

M,  =  AH(e)RliA(0)  (2.4a) 


M2  =  Am(0)R;J  A(8)  (2.4b) 


Note  that  in  accordance  with  the  minimax  principle 
underlying  Semi-coherent  MVDR  MFP  espoused  previ¬ 
ously,  5  is  taken  to  be  that  value  which  maximizes 

^SemiC^  >5).  Since  ^Semi ( r  5)/g„  = 

(S+N)/N~='S/N  +  1,  this  is  equivalent  to  maximizing  the 
beamformer  output  SNR.  Defining  P  =  M7'6,  we  may 
alternatively  define  Ssoni(7)  as 


Ssotu(?)  =  Maximum 
P 


0HAH(6)R7jA(0)P 

PhAh(0)R7?A(8)P 


(2.5) 


The  maximizing  p  is  the  generalized  eigenvector  associ¬ 
ated  with  the  largest  generalized  eigenvalue  of  the  LxL 
matrix  pencil  {AH(0)R7^A(0) ,  Ah(0)R7i2A(0)).  In  turn, 
Ssani(^  the  largest  generalized  eigenvalue  of 
(Ah(0)R7«A(0)  .  Ah(0)R-u2A(0)). 


3.  Comparison  of  Semi-coherent  MVDR  MFP 
With  Full-coherent  MVDR  MFP 

Let  c(7)  denote  the  Lxl  vector  composed  of  the 
relative  complex  amplitudes  of  the  L  multipath  arrivals 
for  a  candidate  source  location  7;  c(7)  is  normalized 
such  that  c”(7)c(7)=  1.  The  Full-coherent  MVDR 
MFP  beamformer  for  a  candidate  search  location  7  is  the 
solution  to  the  single  constraint  MVDR  problem 

Minimize  E{ |w*(7)x(n)|2)  =  wH(7)R„w(7) 

w<7) 


w(7)  = 


_ 1 _ 

cH(7)AH(0)R7iA(0)c(7) 


R*»  A(0)c(7)(3.2) 


Similar  to  Semi-coherent  MVDR  MFP,  we  define  the 
value  of  the  Full-coherent  MVDR  MFP  ambiguity  sur¬ 
face  at  a  point 7,  denoted  Sfuu(7),  as  wH(7)Rltw(7)  / 
wH(7)w(7)  with  w(7)  given  by  (3.2).  Again,  this  is 
proportional  to  the  output  (S+N)/N  with  the  set  of 
weights  described  by  (3.2).  Substituting  (3.2),  we  obtain, 
after  some  algebraic  manipulation. 


Sjh,ii(?>  = 


cH(7)AH(8)R7lA(e)c(7) 

ch(7)ah(0)R72a(0)c(7) 


(3.3) 


Note  that  the  expression  for  Spui^r)  is  very  similar  to 
the  expression  for  the  objective  function  in  the  defining 
expression  for  Sscmi(^)  in  (2.5). 

We  now  examine  the  performance  of  Semi- 
coherent  MVDR  MFP  in  the  case  of  a  single  source  at7s 
and  no  imerferers.  We  will  assume  the  L  multipath 
arrivals  to  be  100%  correlated.  In  this  case,  the  covari¬ 
ance  matrix  has  the  asymptotic  form 

R„  =  aiA(0s)c(7s)cH(7$)AH(0s)  +  o2I  (3.4) 

where  a2  is  the  sum  of  the  square-amplitude  of  each  of 
the  L  multipath  arrivals.  We  want  to  examine  the  value 
of  Ssemi(7)  at7=7s  when  R„  has  the  asymptotic  form 
in  (3.4).  In  our  analysis  we  will  need  expressions  for  R^I 
and  R7?.  It  is  easy  to  verify  that 

R«  =  y--J-A(0  )c(7s)ch(7s)Ah(0  )  +  — Vp1  (3.5) 

A-l  CXs  -s  ~ *  G„ 


where 


Oj  =  cH(7I)AH(8s)A(es)c(7s)  =  1|A(Qs)c(7s)|!2  (3.6a) 

X1=a^as+o2  (3.6b) 

and  is  the  projection  operator  onto  the  orthogonal 

complement  of  the  1-D  space  spanned  by  A(0  )c(7s) 
such  that 

P1A(es)c(7s)  =  0  (3.7a) 

ch(7s)Ah(0  )PiA(0  )c(7s)  =  0  (3.7b) 

—5  — a 

Squaring  the  expression  in  (3.5),  we  obtain,  after  some 
algebraic  manipulation,  a  similar  expression  for  R^: 

RI?  =  A(0  )c(7s)cH(7s)AH(0  )  +  -4-P1  (3.8) 

M  «S  -*  -1  O. 

Sstmi(7J )  is  the  largest  generalized  eigenvalue  of 
(AH(0j)Ri[  A^)  ,  Ah(0j)R^A(05)).  To  evaluate  this, 

we  first  observe  (3.9) 

lAH(®j)R»I A^))  c(7s)  =  \i  (AH(0J)R« A^)}  c(7s) 

where  we  have  employed  the  asymptotic  forms  of  R„ 
and  in  (3.5)  and  (3.8),  respectively,  and  invoked  the 
property  in  (3.7).  (3.9)  dictates  that  c(7s)  is  a  general¬ 
ized  eigenvector  of  {AH(0  )R^A(0  )  , 

Ah(0  )R;2A(0  ))  associated  with  the  generalized  eigen- 

value  where  A-!  is  defined  in  (3.6b).  We  now  show 
that  the  remaining  L-l  generalized  eigenvalues  are  equal 
to  the  noise  power  cr2. 

Let  ej,  i=l,...JL,  denote  the  generalized  eigenvec¬ 
tors  of  {AH(0s)RliA(0s)  ,  Ah(0s)R"2A(0s)).  Since 
AH(9  )RliA(0  )  and  AH(0  )R;2A(0  )  are  both  Hermi- 
tian,  the  generalized  eigenvectors  are  orthogonal  in  the 
metric  AH(0  )R;2A(0  ).  i.  e..  e"AH(0  )RI,2A(0  )e.  =  0 

for  i*j.  Since  c(7,)  is  a  generalized  eigenvector,  it  fol¬ 
lows  that 

cH(7,)AH(0j)RI2A(0s)ej=O  j= 2 . L  (3.10) 

Substituting  the  expression  for  R^J  in  (3.8)  above  and 
exploiting  the  property  in  (3.7)  yields  the  following 
simplification  of  (3.10): 

cM(i;)AM(0|)A(«,)ei  =  O  j=2 . L  (3.11) 

Invoking  the  asymptotic  -forms  of  R7i  and  Rif  m  (3.5) 
and  (3.8),  respectively,  and  the  properties  in  (3.6)  and 
(3.11),  we  find,  after  some  manipulation,  that  for 
j=2,...,L, 


{A,,(«i)R7iA(0s))  eJ  =  o^{AH(es)R7jA(9j))  e,  (3.12) 

Thus,  since  a\  otj  +  cr2  >  a2,  the  largest  generalized 
eigenvalue  of  (AM(0j)R7i A^)  ,  AH(0i)R72A(Gi)),  and 

hence  tire  value  of  Ssemi(*i)>  is  X]  =  o2  dj  +  a2,  where 
otj  is  defined  in  (3.6a).  In  addition,  the  corresponding 
generalized  eigenvector,  and  hence  the  value  of  p  which 
maximizes  the  objective  function  in  (2.5),  is  c(7s). 

It  follows  from  previous  observations  that  the  sig¬ 
nal  plus  noise  to  noise  ratio,  (S+N)/N,  achieved  with  both 
Semi-coherent  MVDR  MFP  and  Full-coherent  MVDR 
MFP  at  the  point  7S  is  Sscmi(7j  )/  or2  and  SFuU(7s  )/  o2. 
respectively,  and  that  these  are  equal  if  p  =  c(7s)  in  (2.5). 
The  previous  development  demonstrated  that  in  the 
asymptotic  case  p  =  c(7s)  in  the  case  of  a  single  source, 
assuming  "adequate"  SNR,  and  that  Ss«ni(7,)  = 
o2  otj  +  o2.  Since  (S+N)/N  =  S/N  +  1,  it  follows  that  the 
asymptotic  SNR  achieved  with  both  Semi-coherent 
MVDR  MFP  and  Full-coherent  MVDR  MFP  at  the  point 
7S  is  (OjCts  +  a2)/a2  -  1  such  that 

2 

SNRFu11(7s  )  =  SNRScmi(7s  )=  -^-||A(0  )c(7s)||2  (3.13) 

where  ||x||  denotes  the  2-norm  of  the  vector  x.  Again,  the 
fact  that  Semi-coherent  MVDR  MFP  yielded  the  same 
output  SNR  at  the  true  source  location  as  that  achieved 
with  Full-coherent  MVDR  MFP  hinged  on  p  being  equal 
to  the  vector  of  normalized  complex  amplitudes  c(7s ). 
We  only  proved  this  to  be  true  in  the  asymptotic  case 
with  a  single  source.  In  general,  p  will  only  be  approxi¬ 
mately  equal  to  c(7s  ).  This  is  the  reason  we  refer  to  the 
method  as  "Semi-cohereni"  in  contrast  to  Full-coherent. 

4.  Incoherent  MVDR  MFP 

Krolik,  Lynch,  and  Swingler  [1]  also  take  a  multi¬ 
rank  signal  approach  to  the  source  localization  problem. 
For  each  point  on  the  search  grid,  they  computing  a 
multiply-constrained  MVDR  bcamforming  weight  vector 
according  to  (2.2)  with  5  equal  to  I,  a  vector  composed 
of  all  ones.  The  value  of  the  "Incoherent"  MVDR  MFP 
ambiguity  surface  at  "7  may  be  computed  according  to 
(2.3)  with  8  =  1,  i.  e..  51^(7)  =  SSoni(7  ;  1). 

5.  Simulations 

A  very  simple  ocean  model  was  employed  in  the 
simulations  to  illustrate  fundamental  aspects  of  each  of 
the  three  versions  of  MVDR  MFP  without  getting  lost  in 
modeling  issues.  In  this  model,  an  iso-velocity  ocean  is 
assumed  with  a  sound  speed  of  1500  m/s  and  a  nominal 


depth  of  4500  m.  The  receiving  array  is  vertically 
oriented  and  composed  of  30  sensors  equi-spaccd  by  a 
half-wavelength  (7.5  m)  at  the  source  frequency,  100  Hz; 
the  array  center  is  situated  at  a  depth  of  100  m.  A  single 
source  is  located  at  a  range  of  20  Km  and  a  depth  of  1 200 
m;  it  is  assumed  that  there  are  L=  1 1  dominant  ray  paths 
between  source  and  array.  The  1 1  multipath  signals  were 
pairwise  95%  correlated  and  each  was  attenuated  in 
accordance  with  three  factors:  bottom  lasses,  auenuation 
proportional  to  path  length,  and  cylindrical  spreading  fac¬ 
tor  [2],  The  amplitude  of  each  multipath  arrival  was 
selected  so  that  the  sum  of  the  square-amplitude  of  each 
multipath  arrival  is  equal  to  what  it  would  be  if  the  mul¬ 
tipath  signals  arrived  with  equal  strength  having  an  SNR 
of  0  dB  per  element  The  wavefront  associated  with  each 
multipath  arrival  is  modeled  as  being  planar  across  the 
face  of  the  array.  Finally,  the  noise  was  modeled  as  being 
independent  from  sensor  to  sensor  and  of  equal  power. 
This  simulation  scenario  is  similar  to  that  employed  by 
Krolik,  Lynch,  and  Swingler  in  [1]  except  that  in  their 
case  the  array  was  horizontally  oriented  and  the  mul¬ 
tipath  arrivals  were  of  equal  strength  with  0  dB  SNR  per 
element. 

Figure  1  compares  the  performance  of  the  Full- 
coherent,  "Incoherent",  and  Semi-coherent  versions  of 
MV  DR  MFP  in  the  simulated  ocean,  source,  and  array 
environment  described  above.  (Thus,  any  simplification 
or  deficiency  in  the  simulation  scenario  was  common  to 
the  data  supplied  to  all  three  algorithms.)  Note  that  peak 
heights  are  in  units  of  dB  and  are  relative  to  the  lowest 
point  on  the  surface  being  0  dB.  For  both  the  no 
mismatch  and  mismatch  cases,  the  same  sample  covari¬ 
ance  matrix  formed  from  250  snapshots  was  used  for 
each  of  the  three  versions  of  MVDR  MFP.  In  the 
mismatch  case  each  version  of  MVDR  MFP  operated 
under  the  assumption  that  the  ocean  depth  was  4500  m 
when  it  was  in  fact  4498  m. 

Three  main  conclusions  may  be  drawn  from  the 
simulation  results:  First.  Full-coherent  MVDR  MFP  is 
extremely  sensitive  to  mismatch.  A  2  m  error  in  the 
assumed  ocean  depth  caused  roughly  a  IS  dB  drop  in  the 
peak  of  the  Full<oherent  MVDR  MFP  ambiguity  sur¬ 
face.  Second,  the  performance  of  ’Incoherent"  MVDR 
MFP  is  substantially  degraded  relative  to  that  of  Full- 
coherent  MVDR  MFP  in  the  no  mismatch  case.  The  peak 
of  the  "Incoherent"  MVDR  MFP  ambiguity  surface  with 
no  mismatch  is  roughly  15  dB  less  than  that  of  the 
corresponding  Full-coherent  MVDR  MFP  ambiguity  sur¬ 
face.  Third,  the  performance  of  Semi-coher?nt  MVDR 
MFP  is  both  relatively  robust  to  mismatch,  with  respect 
to  error  in  the  assumed  ocean  depth,  and  comparable  in 
performance  to  Full-coherent  MVDR  MFP  under  no 


mismatch  conditions. 

These  observations  may  be  explained  by  examin¬ 
ing  the  extent  of  changes  in  the  signal  arrival  parameters 
when  the  ocean  floor  is  raised  by  2  m.  Note  that  this 
causes  a  negligible  change  in  the  arrival  angles  of  and  the 
relative  attenuations  amongst  the  1 1  multipath  signals.  In 
contrast,  since  2  m  is  a  significant  fraction  of  the 
wavelength,  15  m,  the  relative  phases  amongst  the  11 
multipath  arrivals  change  dramatically.  This  extent  to 
which  the  phase  of  a  given  path  changes  depends  on  the 
number  of  bottom  (B)  bounces  (and/or  top  (T)  bounces) 
incurred  en  route  from  source  to  array.  For  example,  the 
2  m  change  in  ocean  depth  causes  a  34.5°  change  in  the 
phase  of  the  single  bouom  bounce  path  B  while  causing  a 
227.2°  change  in  the  phase  of  the  triple  bottom  bounce 
path  BTBTBT. 

Full-coherent  MVDR  MFP  assumes  complete  and 
perfect  knowledge  of  the  emitter  steering  vector,  which  is 
a  linear  combination  of  the  path  steering  vectors.  The 
path  steering  vectors  are  relatively  unchanged  when  the 
ocean  floor  is  raised  by  2  m.  However,  the  path  combin¬ 
ing  coefficients,  which  depend  on  the  relative  phases 
amongst  the  multipath  arrivals,  change  dramatically 
causing  a  catastrophic  loss  in  the  performance  of  Full- 
coherent  MVDR  MFP. 

Interestingly,  comparing  Figures  1(c)  and  1(d),  it  is 
observed  that  "Incoherent"  MVDR  MFP  yielded 
enhanced  performance  in  the  mismatched  case  when  the 
ocean  floor  was  raised  by  2  m.  Although  this  seems 
counter-intuitive  at  first,  note  that  in  the  single  source 
case  the  signal -only  (noise-free)  output  of  the  beam 
formed  with  the  "Incoherent"  MVDR  weight  vector  is  the 
direct  sum  of  each  multipath  signal  as  measured  at  the 
array  center.  Fortuitously  for  "Incoherent"  MVDR  MFP, 
it  turns  out  that  for  these  parameters,  this  vector  sum  was 
larger  when  the  ocean  floor  was  at  4498  m  then  at  4500 
m.  In  a  sense,  with  the  ocean  floor  at  4498  m  the  mul¬ 
tipath  signals  were  more  "in-phase"  than  they  are  with 
the  ocean  floor  at  4500  m.  This  demonstrates  the  strong 
dependence  of  "Incoherent"  MVDR  MFP  on  the  relative 
phases  amongst  the  multipath  arrivals. 
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Figure  1.  Ambiguity  surfaces  generated  via  the  Full-Coherent,  'Incoherent"  and  Semi-coherent  versions  of  MVDR  MFP 


